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Abstract. Can we reduce Quantum Field Theory (QFT) to a quantum computation? Can physics 
be simulated by a quantum computer? Do we believe that a quantum field is ultimately made 
of a numerable set of quantum systems that are unitarily interacting? A positive answer to these 
questions corresponds to substituting QFT with a theory of quantum cellular automata (QCA), and 
the present work is examining this hypothesis. These investigations are part of a large research 
program on a quantum-digitalization of physics, with Quantum Theory as a special theory of 
information, and Physics as emergent from the same quantum-information processing. A QCA- 
based QFT has tremendous potential advantages compared to QFT, being quantum ab-initio and free 
from the problems plaguing QFT due to the continuum hypothesis. Here I will show how dynamics 
emerges from the quantum processing, how the QCA can reproduce the Dirac-field phenomenology 
at large scales, and the kind of departures from QFT that that should be expected at a Planck- 
scale discreteness. I will introduce the notions of linear field quantum automaton and local-matrix 
quantum automaton, in terms of which I will provide the solution to the Feynman's problem about 
the possibility of simulating a Fermi field with a quantum computer. 
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1. INTRODUCTION 

Quantum Theory (QT) is a theory of information: this has been shown in Ref. [1], where 
QT is derived from six fundamental assumptions about how information is processed. 
Ref. [1] achieves the goal of a long-term program initiated in 2003, which has been the 
topic of all my talks at the Vaxjo conferences on quantum foundations, since my first 
proposed set of axioms for QT presented in 2005 [3] (see Refs. [4, 5] for an historical 
excursus). The program is certainly not finished, and one of the next tasks is to find other 
equivalent sets of axioms. One fact, however, is definitely established: QT is a theory of 
information, and whatever set of axioms we consider, they must be equivalent to a set of 
rules for information processing. 

The new informational point of view has far-reaching consequences on our way of 
understanding physics, and has already shown its full power in greatly simplifying 
numerous derivations of theoretical results of quantum information theory, reducing 
some theorems to tautologies. As for any new really fundamental idea, the possible 
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logical consequences and the new reasonings based on the new principles are virtually 
unbounded, and many new interesting results are expected. 

The six principles for QT have nothing of "mechanical" nature, since the QT de- 
rived in Ref. [1] is just the general abstract theory of systems — Hilbert spaces, alge- 
bra of observables, states, operations, and effects — and has no bearing on the theory of 
mechanics — particles, dynamics, and quantization rules: this is why I deliberately used 
the words Quantum Theory in contrast to Quantum Mechanics. However, Quantum Me- 
chanics is only a facet of the many sided Quantum Field Theory (QFT) — precisely its 
restriction to a fixed number of particles — and QFT, strictly speaking, is itself a theory 
of systems. It can be objected that also QFT needs quantization rules, but many would 
agree that we should avoid them in a autonomous foundation of QFT: on the contrary, it 
is classical mechanics the theory that should be derived from QFT with a classicalization 
procedure. 

But, how can we formulate a QFT autonomous from classical theory? The first 
thing to do is to take QT as "The Theory" ab initio. The motivation is simple: QT 
as a theory of systems is solid and autonomous. The idea is then to take QFT as the 
theory of infinitely many quantum systems. Somebody may object that at least space 
and time are not systems, and would be missing in such a theory. And what about 
Relativity? The answer to this objection is that space and time simply must be taken 
as non fundamental, but as derived notions: metric should emerge from the purely 
discrete geometry of a network of quantum systems in interaction. With metric emerging 
from pure topology, the only two additional principles will be: topological homogeneity 
and locality — the informational equivalent of the universality of the physical law. The 
Minkowski metric must follow only from locality and causality, the latter being the first 
of the six informational principles for QT [1]. 

In the above scenario QFT becomes a theory of the quantum cellular automata (QCA) 
that are local, translational invariant, and reversible. This digitalization of QFT is the 
natural settings for the holographic principle, with the QCA-based theory becoming a 
more fundamental layer for QFT, which in turn should be regarded as only an effective 
theory, and we will then interpret the space-period of the automaton a to be the Planck 
length. One of the good news about a QCA-based QFT is that, being naturally a lattice 
theory, it will also have the great bonus of avoiding all problems due to the continuum 
which plague QFT (ultraviolet divergences, the Feynman path integral, non localizability 
of measurements, and many more). The bad news is that the digital theory will likely 
miss some of the simplicity of the continuum, and we must seek simple rules for 
interpolating the digital with the analog descriptions, in order to make easy predictions 
of first-order corrections to QFT and breakdown of covariance at the Planck scale, and 
being able to make comparison e. g. with doubly special relativity [6]. 

One striking feature of the computational paradigm is that large-scale Lorentz co- 
variance is a free bonus. As a matter of fact, Lorentz covariance must emerge from the 
computation if the QCA is capable of simulating physics. And, as shown in Ref. [7], the 
Dirac equation (at least for space dimension D = 1) is just the equation describing the 
free flow of quantum information. In Ref. [5] I showed how space-time can emerge from 
a topologically homogeneous causal network (as a quantum circuit) in (1+1) dimensions, 
by building up foliations over the network, the Lorentz time-dilation corresponding to 
an increased number of leaves within the tic-tac of an Einsteinian clock, and space- 



contraction resulting from the inversely decreased density of events per leaf. In Ref. [8] 
with A. Tosini we have shown how a digital coordinate system can be obtained by the 
analogous of the Einstein's procedure over the causal network, with signals sent back- 
and-forth to events from an observer clock. This work is concerned with pure causality, 
without further specification of the nature of the causal relation. For space dimension 
D > 1 one has the Weyl tiling problem [9], stating the impossibility of emergence of an 
isotropic metrical space from a discrete geometry. 2 In Ref. [1 1] I proposed that the same 
quantum nature of the causality relation provides a way out, with superposition of paths 
restoring the isotropy of maximal speed of propagation of information, and I will discuss 
this here shortly in the concluding section. 

In Refs. [7, 12, 13] I introduced the QCA that implements the Dirac equation in (1+1) 
dimension. There are already simulations of the Dirac field in the literature, noticeably 
the lattice-gas realization of Ref. [14], however, they are are not QCA implementations, 
and work only in 1st quantization and are designed for taking the continuum limit. In 
the same Refs. [7, 12, 13] I showed that the euristical derivation of the Dirac equation 
as free flow of information along the network also provides a kinematical definition of 
the inertial mass in terms of the slowing down of the information propagation due to 
the zig-zag motion over the causal network, whereas the Planck constant becomes just a 
conversion factor between the informational mass and the customary mass. In the QCA 
realization of the Dirac equation unitarity in conjunction with discreteness introduces 
a renormalization of the speed of light, corresponding to a refraction index of vacuum 
that is monotonically decreasing vs mass, with complete halt of the information flow at 
the Planck mass — a new phenomenology coming from the informational point of view. 
I also showed how Fermi fields can be eliminated, having a QCA working locally on 
qubits: for D = 1 space dimension this solves the problem raised by Feynman [15], who 
questioned the possibility of simulating a Fermi field by a quantum computer with local 
interactions. 3 Here in Sect. 3.3, after reviewing this solution for D = 1, 1 will also provide 
the solution for D = 2, which also generalizes to D = 3. 

We have now a fully ab-initio quantum field theory made with a QCA. The classical 
field Hamiltonian is derived from the QCA as emergent, and this also opens the route to 
the mentioned classicalization procedure. In order to address the Feynman problem in 
precise terms I will also introduce the convenient notions of linear quantum field cellular 
automaton (LQCA) and of local-matrix quantum cellular automaton (MQCA). I will 
prove that the vacuum of the MQCA is unique and factorized into single qubits states 
for D = 1 , whereas for D = 2 additional qubits are needed corresponding to auxiliary 
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Majorana fields, whose state is entangled. I will also provide analytical relations between 
the "emergent" and the "time-interpolating" Hamiltonians, and devote an entire section 
to the Quantum Random Walks of the LQCA, namely its first-quantization, providing 
the first numerical evaluations for one and two particle states in the Planck regime. I will 
conclude the paper with some open problems and a brief discussion about possible ways 
out. 



2. THE FIELD-LINEAR QUANTUM CELLULAR AUTOMATON 

We will consider quantum fields on a one-dimensional lattice Z. We will use the follow- 
ing notations for the field: 

1 . (j> n in G Z) denotes a generic scalar (Boson or Fermion) field; 

2. J3 = {<!>„}, a in a finite set, denotes the vector field whose component generate 
the local algebra of the automaton at n\ 

3. for the specific case of Dirac field the letter is substituted by the letter i//; 

4. % denotes a generic scalar anticommuting field. 

Later on we will also consider fields on a D-dimensional lattice with D > 1 , e. g. lP , and 
will use the boldface notation n for the labeling on the lattice. We consider the special 
case of quantum cellular automaton, whose algebra evolution is assigned by a linear 
evolution of a quantum field on the lattice. 4 We now focus on the case of D = 1, and 
consider D > 1 in Subsect. 3.5 and in the concluding section. 

In a linear quantum field cellular automaton (LQCA), the (formally unitary) operator 
U of the evolution transforms the field as follows 



t(t+t) =U^(t)U = Ut(t), 



(1) 



where t is a multiple of the time t of each step of the automaton. Since the evolution must 
preserve the (anti)commutation relations for the field, the unbounded matrix U := \\Uij\\ 
must be itself unitary. The inverse evolution thus is = U<j>(t+ i )U* = U'0(f+ t). 

A LQCA evolves the operator algebra locally. This corresponds to having only a finite 
number of nonvanishing elements in the rows and columns of the matrix U, namely U is 
a band matrix. An example of LQCA is depicted in Fig. 1 with local algebra generated 
by the Fermi-field vector 
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4 Considering boundary conditions on a bounded lattice does not affect derivations as long as we contem- 
plate evolutions for finite numbers of time-steps of Fock states, namely states localized over a quiescent 
"vacuum" (see the following), and we can take the evolution as formally unitary, even for unbounded 
lattice. 
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FIGURE 1. (Left) Illustration of the linear cellular automaton for simulating the Dirac-Weyl equation in 
one space dimension. (Right) Margolus scheme for the automaton, namely the quantum circuit mimicking 
the Dirac equation. 



As proven in Ref. [16] a cellular automaton can be always decomposed into a finite num- 
ber of rows of independent gates connecting field operators in neighboring locations: this 
is the Margolus scheme [17] (see Fig. 1). The cellular automaton is thus described by 
a quantum computation that is homogeneous both in space and time, corresponding to 
a causal network with square lattice. On such lattice the Feynman path-integral is well 
defined, and corresponds to composing the evolution from many gates in deriving the 
input-output relation of fields. The path-sum rules for the evolution have been given in 
Ref. [7]. 



3. THE DIRAC LQCA 

The (spinless) Dirac equation in 1 space dimension in the Weyl representation is 

d t \if=(-ca 3 d x -io l co)}if 1 



(3) 



where co = cA _1 , X = h/mc is the Compton wavelength, and with vectors | f) and 
1 1) representing right-handed and left-handed fields, respectively, whereas particle and 
antiparticle fields correspond to vectors |+) and |— ), respectively. In matrix form we 
have 
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The field l/A satisfies the anti-commutation relations [\]/ a (x), Y (y)^]+ = &ap8( x ~ y)> 
xj6K denoting values of the space coordinate. We now want to derive a LQCA 
corresponding to the discretized version of the Dirac equation. We label the field \jf n 
with a discrete index nGZ corresponding to discrete positions na, with a the periodicity 
of the LQCA. The LQCA version of the Weyl equation would be the same as Eq. (3) 
with the partial derivatives substituted by finite differences (it is convenient to consider 

finite differences in the symmetric form d := h(d + — <9_), with d±f(n) := f(n ± 1)). In 
the quantum automaton version we will take for the speed of light c the causal speed 
c = a/ 1, namely the maximal speed of propagation in the automaton. Then, due to the 
unitarity, a discretized version of the Dirac partial differential equation (4) is possible 



only if we allow for a renormalization c — > £ c [7] (see also the following), namely 
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where d t 
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a [ d (with the shift d± in the appropriate discrete variable), 



t denoting the execution time of the LQCA. The time-difference in the LQCA (5) 
corresponds to the difference between the unitary matrices U and U^, given by 
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where 



c = cot = a/A, s = £, 
and with unitarity implying the identity 

c 2 + s 2 = 1. 

Using Eq. (7) and parameterizing c and s by an angle 6, one has 

s = sin = £ = 
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which shows that is a mass-dependent vacuum refraction index which is strictly 
greater than 1 for nonzero mass, monotonically increasing versus m, and becoming 
infinite at m = m. For a the Planck length m is the Planck mass, and the automaton 
becomes stationary (i. e. there is no propagation of information) at the Planck mass: this 
interesting violation of dispersion relation has been presented in Ref. [7], and is due only 
to discreteness in conjunction with the unitarity of the automaton (see Fig. 2). 



3.1. Margolus scheme for the Dirac LQCA 

The Dirac automaton (5) can be achieved with the Margolus scheme in Fig. 1. 
The form of the gates A and B can be derived assuming, without loss of generality, 




FIGURE 2. The inverse vacuum refraction index £ versus the mass m of the Dirac field. The mass scale 
is given by the Planck mass m = a denoting the period of the automaton (from Ref. [7]). 



unit determinants |A| = |B| 

A21 =Mi = Tis, and A 22 = 



= 1. Then, one has [7] B\\ — B22 = 0, B12 = #21 = ±*\ 
An = q=c. In the following, we will adopt the solution 
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3.2. Emergent Hamiltonian for the LQCA 

The LQCA has no Hamiltonian: all unitary transformations are far from the identity. 
The Hamiltonian, becomes an "emergent" notion: it can be written in terms of the LQCA 
unitary matrix as follows (we remind the vector notation in Eq. (2)) 

H = ^(V-^)\if. (11) 

It is easy to show that one has 

ihd t \ir=[\if,H] (12) 

and in this sense H is an Hamiltonian associated to the LQCA. It is formally iden- 
tical to the classical field Hamiltonian, which gives the field equation via Poisson 
brackets. Notice that Eq. (11) would also hold using the non Hermitian Hamilto- 
nian H = iht~ l y/^ (\J\ — V^Xjf satisfying Eq. (12) for a halved-time finite-difference 
derivative, with Ui and U2 the unitary transformations associated to the two rows of 
gates in the Margolus scheme. Eq. (12) provides a three-point automaton evolution- 
rule y{t + t) = y(t - t) - 2i ltr l [^(t),H], which can be time-reversed as y(t - t) = 
yf(t+ t) + 2i ifC x [\j/(t),H]. Thus the automaton invertibility is due to the existence of a 
three-point updating rule, and not to Hermiticity of the Hamiltonian: the fact that Hamil- 
tonian can be chosen Hermitian is a consequence only of time-homogeneity of evolution 
(the association of reversibility of a cellular automaton with a three-time updating has 
been first noticed in Ref. [18]). 

The mapping between Hamiltonian and unitary evolution operator for the LQCA can 
be considered the discrete version of the operator exponential mapping. One has 

U = exp[-/sin- 1 (Ht/^)] (13) 

in terms of the matrix 

H:=-(U-U t ) = [[ V A,//], V A t ] + , (14) 

corresponding to 

// = V f HyA. (15) 

This Hamiltonian is different from the time-interpolating Hamiltonian H defined 
through the identity 

U =:exp(-z//t//?), (16) 



which instead is given by 



(17) 



and is generally highly non local. 

3.3. The Feynman problem: simulating Fermi fields with a quantum 



The algebra generated by the field in a given point n of a lattice, e. g. n G Z , is 
what is usually referred as "local" or "quasi-local" algebra in the literature, and this is 
what we called local algebra in the definition of the field quantum cellular automaton. 
For Fermi fields, however, this locality notion is artificial, since the field anticommutes at 
different positions and is itself non observable, as opposed to the case of a Bose field. As 
mentioned in the introduction, the possibility of simulating a Fermi field by a quantum 
computer — shortly qubit-izing Fermi fields — is still an open a problem, and was raised 
by Feynman thirty years ago [15]. The solution to this problem is very relevant, since 
it would respond also to the general question (also raised by Feynman in the same Ref. 
[15]) whether Physics can be simulated by a universal (quantum) computer, this is also 
the statement of the Deutsch's version of the Church-Turing principle [19]. 

Let's state the Feynman problem more precisely. Since the case of Bose field is 
obvious (commutation is trivially achieved by tensor product of local algebras), I state 
the problem only for the Fermi case. We will call a quantum cellular automaton with 
algebra given by the tensor product of local algebras that are finite-dimensional (with 
dimension independent on the size of the automaton in the bounded case) local-matrix 
QCA (MQCA). We remind that the locality of the evolution is a crucial part of the 
definition of QCA. With the above in mind, we now raise the fundamental problem 

Feynman problem: Is it possible to find a local-matrix QCA for any Fermi LQCA? 

For LQCA in D = 1 space dimension a positive answer has been given in Ref. [7], which 
will be shortly reviewed here. Thereafter I will also provide a solution for D = 2 case, 
which can be simply extended also to D = 3. 

The Fermi field can be put into one-to-one correspondence with an algebra generated 
by a tensor product of local Pauli algebras by a relation of the form 



where the "phase factor" 4>(n) ^ is a nonlocal functional of <p n . Analytical expressions 
for 4>(n) has been given by Jordan and Wigner in Ref. [20] for D = 1, and for D = 2 
and D = 3 by several authors in Refs. [22, 23] and Ref. [24], respectively. Unfortunately, 
as we will see in the following, these expressions cannot be used to solve the Feynman 
problem for D > 1: to such purpose, a selfadjoint 4>(n) is needed, and this new JW 
relation will be presented here. 
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3.4. Solution of the Feynman problem for space dimension 1 



Let me briefly review the solution of Ref. [7]. The Jordan and Wigner [20] form of 
the phase factor 4>(n) is given by 



®(n) =exp ( iTtY,<Pk<Pk 

k<n 



and identity (18) is inverted as follows 



<pj = n.(-°k)°j ■ 

k<j 



(19) 



(20) 



We see that the Pauli-local algebras are non-local in the field-local algebras, and vicev- 
ersa the field-local algebras are non-local in the Pauli-local algebras. Using the following 
identity that holds for both commuting and anticommuting z, z' 
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it is easy to check that 
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with 

A n = exp 



t,„+ 



B n = exp 



(23) 

Using the JW identity (19) with = (p2 n and \\f n = (fen+i, the unitary operators of the 
gates rewrite as follows 

A n = exp [-W {C^ n _ x ctn + °2n-l°2n) ] . B n = [f + °2n°2n+\) ] ■ 

(24) 

From the above construction it is also clear that any LQCA in (1 + 1) dimension has 
a MQCA realization, since it has a Margolus realization with two rows, and each 
gate unitary operator can be expressed as exponential of bilinear form of the field in 
neighboring positions, say n and n + k, which can be written as tensor products 
of local matrices in all positions starting from n and ending at n + k. Indeed, the JW 
construction leads to the "string" identity 

<p]+i<Pj = (-) l orS l jO+ +l , S l j:= J] °l ("string" identity), (25) 

)<k<j+l 

and only for neighboring locations (/ = 1) one has Pauli operators only at the same 
locations of the fields, otherwise one has a string of additional o 3 operators joining the 
two locations. Can we reduce the number of qubits to those only at the same locations 



of the field? An idea for answering to this question is provided by Ref. [28], and 
consists in introducing additional Majorana fields <pj := <pt 4- cpj and <pj := i(<pt — <pj) 
anticommuting with the original Dirac fields, and writing both Dirac and Majorana fields 
in terms of a single JW string of Pauli operators, pertaining Dirac and Majorana fields in 
alternate positions (in our case taking the even locations for the Dirac and the odd ones 
for the Majorana). Then, one has the identity 

^j + 2^jP}+2i + i,j + i = -ioJof +2l aj +l a) +2l+l , pJj = P tj := i<p} <p] , pf j = L (26) 

The idea now is that if all involved in the interaction terms <pj +2 / < Pi commute > 
then, upon jointly diagonalizing them on a common vacuum state, one would have 
<Pj +2 i<Pj = ^ z ^7 (7 j+2/ <7 J+i (7 j'+2/+i' an d the MQCA on that vacuum would be local. 
However, it is easy to check that [Pij^Pff] = [Pij,Pji] = for i 7^ i' and j 7^ /, but 
[Pij,Pn] + = [Pn,Pji] + = 0. Therefore, the construction works if there are interaction 
terms only at a single fixed distance Z, not when there are more distances / of interactions. 

3.4.1. The Vacuum for the Dirac automaton 

There are many definitions of vacuum in QFT, and they are no longer equivalent in 
the case of the quantum automata. In order to build up a Fock space for the LQCA, we 
need to require that the vacuum \Q) be annihilated by the field operator <p n , namely 

(p n \a)=0. (27) 

Notice that in the present context such definition has no apparent relation to the vacuum 
definition as the lowest-eigenvalue eigenstate of the Hamiltonian, even for the emergent 
Hamiltonian (11), since this is not jointly diagonal with U. On the other hand, as we will 
see in Sect. 4, |£2) is left invariant by the unitary evolution. 

In the JW representation the only possible vacuum for the D = 1 Dirac MQCA is 

\a) = ®k\l)k- (28) 

Proof: Write the vector |Q) as =: | l) n \&i) + | t)n\&l), with and (gen- 
erally are unnormalized and non orthogonal) states for the algebra of all qubits at 
location m ^ n. Being 4>(«) an exponential, it cannot annihilate any vector, whence 

= I \.} n ^>(n)\Q.l}, which is null if and only if = 0. Upon repeating the same 
reasoning for all n, the statement (28) follows. ■ 
Notice the useful identity: (p} t \^l) = 0"+|£2). 

3.4.2. Relation with one-dimensional models in statistical mechanics 

It interesting to notice that the Hamiltonian (11) in the JW realization of D = 1 
Dirac field corresponds to a spin model of the form H = — i c £ a -1 o^a^_ l a~_ 2 — 



w Ln CT 2n°2n-i + n - c - It would be interesting to see how an integrable spin model 
(solved e. g. by the Bethe ansatz and/or Yang-Baxter) is translated into a LQCA, in 
which case the LQCA would provide an easy numerical simulation for the evolution 
for discretized time for first-quantization situations as e. g. spin-waves. Notice, how- 
ever, that for example the conventional Hamiltonian for the Heisenberg XY model in 
transverse field corresponds to a Dirac conventional Hamiltonian containing also a four- 
linear interaction term, indicating that the field-linearity of the automaton may not be 
preserved in the correspondence. 

Another interesting relation with condensed-matter theoretical models is the case 
of Hamiltonians of the form H = — hK^jiCj^Cj + cjcy+i), K a coupling constant, 
describing noninteracting electrons hopping on a 1 -dimensional lattice. There is noting 
relativistic about this: however, if you fill the system up to the Fermi energy Ep and focus 
on electrons with energy £ = £ - < £ f and momentum \k — kp\ <C k F , the right and 
left-moving fields associated with this electron will obey the zero-mass Weyl equation 
in 1 space dimension. Then, a Peirls distortion would open an off-diagonal term in the 
energy spectrum corresponding to a mass-coupling between the two fields (see e. g. Ref. 
[25]). 



3.5. Solution of the Feynman problem for space dimension D=2,3 

In Ref. [28] it has been proposed to use Majorana auxiliary fields to eliminate fields 
in Hamiltonian models for D > 1. The idea is to use the construction mentioned in 
Subsect. 3.4, based on covering the D-dimensional lattice with a 1 -dimensional string. 
Unfortunately, as I have shown, the construction works only if there are interaction terms 
at a single given distance, whence it cannot work in the 1 -dimensional covering, since 
there are always more than two neighboring locations in interaction (the coordination 
number is necessarily greater than 2). In addition, the 1 -dimensional covering raises 
the problem that the bulk realization will depend on boundary terms, and would be 
unsuitable for unbounded lattices. The idea of using Majorana fields, however, remains 
valuable, and will be used here. In constructing the Majorana fields, it has been crucial 
to have the JW "phase factor" as formally selfadjoint $(n) = 4>^(n), so that 2 = /, 
which is needed for the cancellation of strings of Pauli matrices. Unfortunately the form 
proposedd in Refs. [22, 23, 24] for D = 2, 3 do not satisfy this requirement. We need 
then a new JW construction. 

In order to have identity (18) connecting commuting Pauli matrices c n with anticom- 
muting fields (p n , one needs to satisfy the identity 

■ (m) 

cp n ®(m)=<i>(m)( Pn e ia » , (29) 
with the phase satisfying the identities 



c4 n) = oi m) +^ «i n) =0. 



(30) 



In order to satisfy identity (29) one needs to adopt the general form for the phase factor 

4>(n) :=exp (^Y, ( Pm ( PmO& ) ^j , (31) 

whereas for D = 2 identity (30) along with the requirement that 4> is selfadjoint, imply 
the following step-function analytical form for oc^ 

(m) a < >\ a n \ J ' if < arg(k : i) < ^, 

(A ; = 2 ^ m-n:i , Q l7Z k = <^ ) (32) 

(ft, if % ^ arg(k : l) < 2n, 

where arg(k : i) is the angle between the vector k and the vector i directed as the x-axis 
in the plane. It is immediate to check that 4>(n)' f ' = $(n) and <3>(n) 2 = /, and that for 
D = 1 we recover the JW construction, whereas for D = 2 the construction corresponds 
to set the origin of the lattice in n and put an identity in each site of the upper half -plane 
(including the positive x axis and the origin), and put a cr 3 on all other sites. For D = 3 we 
need to make a rotation by n around the axis orthogonal to the plane identified by vectors 
n and m, and in order to keep the transformation of Dirac spinors under rotations we can 
use the isospin form of the phase-factor suggested in Ref. [24], but with the angle again 
as a step-function (32), in order to satisfy the additional requirement $(n)' i ' = 4>(n). In 
the following we will focus for simplicity only on the case D = 2. 

We can now add auxiliary Fermi fields e. g. upon redefining the phase factor as follows 



4>(n) := exp 



(33) 



and introduce corresponding Pauli matrices as follows 

T+ = ^4>(n). (34) 

The auxiliary fields anti-commute with the original Fermi fields. 5 We can now introduce 
the Majorana fields with v = 1,2 corresponding to the Pauli matrices T„ and 
introduce the observables P mn := z^^n- It is easy to see that P„ m = /, and for n ^ 

W and m ^ m' one has [Pnm^n'm'] = [^nm^mn'] = [^m,4n] = 0, [Pnm,^nm']+ = 

[^mn 5 ^Vn]+ = 0- Then one has the identity 

(piWmPmn = ~i<?n <? m ?n?m- ( 35 ) 

If we now orient all interaction links 1 and take only +1 when we have both ±1, then 
^n,n+i for varying n and 1 commute each other (as long as we do not have different 
1 that are parallel). One can then introduce a joint vacuum for both types of fields, 
which jointly diagonalizes all the observables in, n +i- It is easy to see that in the qubit 



5 A. Tosini has shown that it is possible to make the Dirac and Majorana fields both commuting and 
anticommuting. 



representation such vacuum is factorized into separate components for o and x algebras, 
and the derivation of Subsect. 3.4.1 still holds for the a part which is then a tensor 
product of 1 1) whereas the T- vector is entangled. On such vacuum, one has the identity 

<Pn<Pn+l + <P,I +1 <Pn = ±(<Pn<Pn+l + <pl +l <Pn)Pn,n+l = ±/(c7+C7 n+ l - C7+ +1 (7 r ") (36) 

The same vacuum diagonalizing the observables in, n +i does not diagonalize also T^x^ +l , 
and we need to work with a doubled qubits. 



4. FIRST QUANTIZATION: QUANTUM RANDOM WALKS 



In this section we come back to the D = 1 Dirac automaton. The number of up-qubits 
S = Em Gn i s a constant of motion. The observable algebra thus can be decomposed 
as srf = 05 = o=2^> and one can consider the evolution within each single subalgebra 
separately. A basis for the single-particle Fock sector is given by the orthonormal set 



{|V„ a )}a=±,n€Z, \K)-=¥^), 



(37) 



or, in tensor notation 



\<) = \n) 



\a). (38) 

Using the invariance of the vacuum U\Q.) = and the identity U\yf n ) = U y^U^Gl) = 
V4U|n)=U|vO>we get 

tti v Ts 0l , , — ic] . , To ol i > 

U\¥ n )= n o IV„-i>+ _ ir n IV«> + n „ IV»+i>, (39) 



or, in tensor notation, 

U = s(e_<g>| t)(t \+e+®\l)(l I) -z'c/tgxr 1 , 



where e± := |n ± l)(n|. The evolution of single-particle states 



It2n), 
I t2n+l), 



a = + 
a = — 



(40) 



(41) 



is a unitary band-matrix with matrix blocks given by Eq. (39). In Fig. 3 as an example 
we give the computer simulation of a Gaussian single-particle state of the form 

\ x i')=N k -^-(|y+)± +: particle,-: antiparticle. (42) 

In the same Fig. 3 we also see a double-slit particle localized state evolving. It is 
noticeable how the Gaussian packet spreads very little compared to localized states. 
One can see that the following states are invariant 



s sin + al, 
c 



£ = \/l — S 2 COS 2 0. 



(43) 



Notice that such states are not normalizable for unbounded LQCA. For small mass (and 
or small "momentum" sin0 = e ) the invariant states approach left and right-handed 
states, whereas at the Planck mass limit they become particle/antiparticle states. 
In the single-particle sector one can define the position and momentum operators 

P a = -mZ\ ¥ Z)l( ¥ Z\,X a = 2a£\ ¥ «}n( V Z\. (44) 



and check the commutation relations [X a ,pP] = ihl a 8 a p, and [X,P] = ihl\. For the 
invariant states (eigenstates of the momentum), one has P| = ^ sin (j) \ Iff) . 
The general form of a two-particle state is (1//+ = (fan, y/~ = (Pin+i) 



|tf 2 )> = £¥g>9>?|n>, * 



,(2) 



(2) 
ji ' 



(45) 



Clearly the (normalized) anti symmetrized Kronecker product of two single-particle 
vectors is a two-particle state, and it is also easy to prove that its unitary evolution is just 
the antisymmetrized Kronecker product of the unitarily evolved single-particle states. 

Finally, regarding the mentioned classicalization procedure, the derivation of the clas- 
sical field theory corresponding to the quantum cellular automaton proceeds as follows. 
We build up the fields from qubits via the Jordan Wigner transformation, deriving the 
classical Hamiltonian via the procedure in Subsection 3.2. Then, the classical particle 
emerges from Gaussian single-particle states, with position and momentum given by the 
expectations of the operators in Eq. (44). The classical trajectory is the "typical path" 
along the quantum network, namely the path with maximum probability (see e. g. Fig. 
3). " 
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FIGURE 3. (Left) The evolution of a single -particle Gaussian packet of the form in Eq. (42) with xo = 0, 
A = 2, k = 8, for 180 time-steps (i. e. t = 180t) and a total dimension of 128, corresponding to 128 qubits, 
half of them for left and half of them for right particles. The red line is the typical path, corresponding to 
the classical trajectory. The parameter c = cos(0) with 9 = 7r/8 here corresponds to m ~ .92m. (Center) 
Details in 3D of the evolution of the single-particle Gauss-packet limited to dimension 64 and for 20 
time-steps. (Right) A "double-slit" particle state I'P) = |(| t+) + | f~) + | ft n ) + | t- n )) for n = 10, 80 
time-steps, and — n/10. 



5. CONCLUSIONS AND OPEN PROBLEMS 



We have explored the quantum cellular automaton introduced in Ref. [7] to implement 
the Dirac equation in (1+1) dimension as a quantum computation. The automaton can be 
regarded as a quantum ab-initio field theory, i. e. made with QT only with the addition of 
the requirements of interaction-locality, translational invariance, and reversibility. From 
the first computer simulations we have seen that at large scales the QCA mimics the 
corresponding QFT theory, whereas interesting departures arise at the Planck scale, for 
high mass and momentum. As mentioned, one of the most urgent open problems is now 
to devise an asymptotic method to evaluate the leading corrections for large number of 
time-steps and states smoothly varying over large number of sites, in order to derive the 
phenomenological modifications occurring at the Planck scale, of the kind of e. g. of Ref. 
[6]. This task is made difficult by the simple fact that discrete automata generally exhibit 
chaos even in the quantum case (see e. g. Ref. [26]): here, however, a simplification is 
provided by the field-linearity of the automaton. 

We have seen how the Feynman problem of simulating the field automaton by a 
local-matrix quantum automaton is solved for D = 2, and similarly can be solved for 
D = 3. Additional problems, however, must be addressed for space dimension D > 1. 
For example, we still must write down explicitly a Dirac automaton for D > 1 (notice 
that generally there exist many automata mimicking the same linear field equation). 
An instance of such automaton is the one provided by Bialynicki-Birula [27], which 
is proved to approach the Dirac equation in the continuum limit a — > (via Trotter's 
formula), but it is unknown how it would behave in the large-scale limit of infinitely 
many steps and smooth delocalized states. A more interesting possibility is to consider 
the emergent Hamiltonian for the automaton version of the QFT theory, and write the 
Hamiltonian as a difference of the unitary matrix of the automaton and its adjoint, 
as in Eq. (12). The Weyl tiling issue for D > 1, on the other hand, clearly pose no 
problem for the automaton realization, since it would persist in the continuum limit, 
but is contradicted by the continuum limit of the Bialynicki-Birula automata. This 
shows that the quantum nature of the causal network plays a crucial role in having 
Minkowski space-time as emergent from the discrete geometry of pure topology. It 
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FIGURE 4. Evolution at varying time-steps N t of a two-particle state made upon anti-symmetrizing the 
two Gaussian packets in collision with xq = ±10, A = 2, k = ±8, 128 qubits. Here m ~ .92m as in Fig. 

3. The plot represents the matrix I'l'iml 2 in Eq. (45), with the vertical and the horizontal axes denoting 
the one-dimensional coordinate of the two particles. The symmetry along the diagonal is the result of 
indistinguishability. 



would be interesting to see explicitly how the quantum superposition of different paths 
on the network restores the metrical isotropy. 
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